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ABSTRACT 

An  analytical  High  Value  Target  (HVT)  acquisition  model  is 
developed  for  a  generic  anti-ship  cruise  missile  system. 
The  target  set  is  represented  as  a  single  HVT  within  a  field 
of  escorts.  The  HVT's  location  is  described  by  a  bivariate 
normal  probability  distribution.  The  escorts  are  repre- 
sented by  a  spatially  homogeneous  Poisson  random  field 
surrounding  the  HVT.  Model  output  consists  of  the  prob- 
ability that  at  least  one  missile  of  a  salvo  acquires  the 
HVT,  conditioned  on  the  number  of  missiles  in  the  salvo 
which  penetrate  the  HVT  area  defense.  In  addition,  the  fall 
of  multiple  penetrators  is  modeled  using  a  conditional 
multinomial  probability  distribution.  The  model's  equations 
are  used  to  solve  for  an  optimal  missile  seeker  range  gate, 
given  a  probability  distribution  describing  the  location  of 
the  HVT  within  the  targeted  formation  at  the  time  the 
missile  commences  its  search.  Included  in  an  appendix  is  a 
time-dependent  model  describing  HVT  location  which  provides 
for  HVT  movement  during  missile  time  of  flight  up  to  the 
moment  of  active  search. 
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I.  INTRODUCTION 

With,  the  advent  of  modern,  sophisticated  anti-ship 
cruise  missiles  (ASCM),  opposing  naval  battle  forces  have  an 
over-the-horizon  (OTH)  strike  potential  which,  when  actual- 
ized by  appropriate  tactics,  will  herald  a  new  era  in  naval 
combat.  As  an  opening  remark  the  above  may  sound  dated,  for 
the  cruise  missile  has  arrived  and  the  technology  is  aboard 
many  combat  fleets  of  the  world.  Nevertheless,  the  implica- 
tion of  next-generation  microelectronics  and  artificial 
intelligence  aboard  the  ASCM  of  tomorrow  presents  the 
forward-looking  military  operations  research  analyst  with 
promising  opportunities  for  tactical  improvement. 
Revolutionary  advances  in  ASCM  technologies  will  provide 
future  systems  with  multi-target  discrimination  ability, 
reduced  soft-kill  susceptibility,  and  other  significant 
improvements.  Hence,  tactical  options  should  be  evaluated 
well  in  advance  of  a  new  missile's  arrival  in  the  Fleet. 

In  an  attempt  to  destroy  an  opposing,  -escorted  high 
value  target  (HVT),  the  decision-maker  possessing  an  ASCM 
OTH  capability  must  choose  when  to  strike,  using  appropriate 
force  ( e. g. ,  total  salvo  size),  a  decision  surely  influenced 
by  targeting  accuracy,  intelligence  regarding  opposing  force 
composition,  and  first-strike  criticality.  Quantitative 
answers  to  the  above  decision  problem  are  by  nature  probabi- 
listic ones.  The  decision-maker  can  be  assisted  by  an  asso- 
ciated estimate  of  HVT  kill  probability  from  a  model. 

The  multi-target  scenario  has  been  modeled  quite  exten- 
sively using  elaborate  Monte-Carlo  computer  simulation  tech- 
niques which,  even  after  considerable  programming  input 
effort  and  program  execution  time,  yield  estimates  that 
often  neglect  the  essential  vagaries  of  the  war-at-sea  envi- 
ronment.   Simulation  has  contributed  much  insight   to  the 


decision  problem  at  hand,  but  is  generally  inappropriate  as 
a  tactical  decision-aid  because  of  its  inflexibility.  The 
so-called  Bernoulli  Trials  model  is  a  much  simpler,  purely 
analytical  probability  model  which  has  been  employed  as  a 
tactical'  decision-aid,  [  Ref .  1]  . 

In  the  Bernoulli  Trials  model,  HVT  defense  is  essen- 
tially limited  by  the  number  of  available  fire-control  chan- 
nels, each  of  which  is  awarded  a  deterministic  or 
probabilistic  number  of  ASCM  kills.  Targeting  information 
is  assumed  to  be  accurate  and  all  ASCM  that  penetrate  the 
HVT  area  defense  have  equal  probability  of  acquiring  any  of 
the  multiple  targets  present.  In  an  extension  to  the  model, 
HVT  acquisition  by  the  ASCMs  is  biased  upward,  thereby 
endowing  the  missiles  with  multi-target  discrimination 
ability.  The  model's  equations  are  solved  iteratively  for 
an  estimate  of  salvo  size  needed  to  inflict  a  firepower-kill 
on  the  escorted  HVT  with  a  specified  probability  of  success. 
The  driving  parameters  are  the  number  of  escorts,  the  effec- 
tiveness of  HVT  defense,  and  the  required  number  of  hits  on 
the  HVT  to  achieve  a  firepower-kill. 

The  Bernoulli  Trials  model  and  its  simulation  counter- 
parts make  simplistic  assumptions  regarding  ASCM  target 
acquisition.  An  accurate  target  acquisition  model  is  an 
essential  part  of  any  model  that  attempts  to  provide  the 
strike  planner  with  a  cost-effective  estimate  of  requisite 
salvo  size,  where  "cost"  refers  to  the  number  of  ASCMs 
employed.  (While  it  is  true  that  the  prudent  tactician  may 
wish  to  launch  "extra"  missiles  as  a  hedge,  it  is  also 
likely  that  missiles  will  be  in  short  supply;  hence,  he 
needs  to  know  what  basic  number  he  is  supplementing  with 
"extras". )  This  paper  offers  a  new,  purely  analytical  HVT 
acquisition  probability  model  designed  to  enhance  existing 
multi-target  ASCM  OTH  models.  It  retains  much  of  the 
simplicity  of  the  Bernoulli  Trials  model,  but: 
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•  Reflects  statistical  fluctuations   in  targeting  accuracy 
in  a  better  way. 

•  Incorporates  the   attacker's   missile   search  area   and 
guides  his  range  gating  decisions. 


II.  THE  MODEL 

A.   DESCRIPTION 

1.   Overview  and  General  Assumptions 

The  ANALYTICAL  HVT  ACQUISITION  MODEL  estimates  the 
probability  that  at  least  one  penetrator  ( an  ASCM  that 
survives  the  HVT  area  defense)  acquires  a  single  HVT  within 
an  escort  field  of  low  value  targets  (LVT).  Additionally, 
the  expected  number  of  penetrator s  that  acquire  the  HVT  and 
her  escorts  ( ie. ,  the  fall  of  shot)  is  estimated.  The 
model's  outputs  serve  as  reasonable  Measures  of 
Effectiveness  of  the  decision-maker's  choices  under  uncer- 
tainty; namely,  given  the  current  tactical  picture  and 
opportunity  to  strike  at  the  enemy,  how  many  missiles  should 
he  employ  and  how  are  they  best  programmed  to  achieve  a 
desired  probability  of  HVT  acquisition? 

The  choices  of  launch  bearing  and  ASCM  search 
programming  are  influenced  by  the  accuracy  and  time  late  of 
the  targeting  information  and  the  missile  time  of  flight  to 
the  target  area.  Where  and  when  opportunities  exist,  anal- 
ysis of  the  offense's  information  state  concerning  the 
targeted  formation  may  be  profitable  even  when  information 
is  scarce  and  uncertain.  The  model  assumes  the  following 
information  is  available  to  the  strike  planner: 

•  A  bivariate  normal  probability  density  function  lor 
other  appropriate  probability  density)  describing  the 
location  of  the  HVT  within  an  array  of  escorts. 

•  An  estimate  of  LVT  density  about  the  HVT. 

•  An  appropriate  set  of  missile  seeker  sensor  Sweep 
Widths. 

•  An  estimate  of  target  acquisition  probability,  given 
target  detection. 
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•  A   simplified  estimate   of  ASCM   reliability  during   the 
search,  detection,  and  acquisition  phases  of  the  attack. 

•  A  velocity-time  random  vector  of  HVT/escort  field  motion 
( optional ) . 

The  missile  firing  strategy  used  in  the  model  is 
that  of  near-simultaneous  firings  ( a  short  firing  interval 
with  respect  to  appreciable  target  movement)  on  a  common 
launch  bearing  calculated  to  intercept  the  mean  of  the  HVT 
location  density.  If  the  model  user  has  a  postulated  motion 
vector  of  the  enemy  formation,  the  HVT  density  can  be 
updated  to  incorporate  the  time  delay  from  the  generation  of 
the  previous  targeting  ellipse  to  the  time  of  arrival  of  the 
salvo  at  the  search  area,  in  which  case  the  salvo  launch 
bearing  should  be  calculated  to  intercept  the  mean  of  the 
updated  HVT  density  at  the  time  of  missile  search. 
Adjusting  the  HVT  density  and  launch  bearing  becomes  neces- 
sary when  missile  time  of  flight  and  enemy  motion  result  in 
significant  displacement  of  the  enemy  formation  from  its 
location  at  the  time  of  launch. 
2.   Modeling  the  ASCM  System 

A  generic  ASCM  system  is  modeled.  Missile  flight  is 
straight  down  the  launch  bearing  through  the  programmed 
search  area.  The  assumption  of  zero  missile  navigation 
error  is  made  since  incorporation  of  the  relatively  small 
errors  associated  with  a  highly  accurate  system  is  not  apt 
to  heavily  influence  model  results.  (Deference  is  made 
however,  to  the  exacting  navigational  accuracy  required  in 
especially  long  range  firings. )  The  missile  seeker  sensing 
device  is  essentially  off  (  ie.  ,  receives  no  information), 
prior  to  ASCM  arrival  at  the  search  area.  An  important 
assumption  concerning  the  offensive  firing  is  the  notion  of 
conditional  independence  between  shots;  each  missile  within 
a  salvo  is  given  the  same  launch  bearing  and  search  program, 
but  retains  functional  independence  ( autonomy)  with  regard 
to  target  acquisition. 
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The  effective  ocean  area  searched  out  by  the  ASCM  is 
a  function  of  operator  input  (range  gate)  and  the  missile 
seeker  sensor  Sweep  Width(s)  given  by  a  Definite  Range  Law 
of  Detection.  The  resulting  area  is  a  rectangle  (or  set  of 
rectangles  when  multiple  Sweep  Widths  are  defined).  In  the 
existing  detection  environment,  individual  targets  may  be 
grouped  by  various  homogeneous  characteristics  (e.g.,  size, 
active  or  passive  signatures,  shape,  motion,  etc. ).  For  the 
purpose  of  detection,  similar  targets  are  subject  to  a 
common  sensor  Sweep  Width.  The  model  may  be  used  in  the 
situation  where  the  HVT  possesses  a  detectable  feature  that 
distinguishes  it  from  the  LVT  class.  For  example,  a 
sensor's  Sweep  Width  for  an  aircraft-carrier  may  be  wider 
than  its  Sweep  Width  for  a  destroyer-escort  because  of  the 
carrier's  relatively  large  size. 

A  target  is  detected  if  it  lies  within  its  associ- 
ated rectangle  of  effective  search  and  is  "swept"  by  the 
missile  seeker  sensor  ( ie. ,  is  in  the  sensor's  field  of  view 
and  within  half  a  Sweep  Width  of  the  search  axis).  The 
model  assumes  the  flightpath  and  the  search  axis  of  each 
missile  are  the  same.  The  model's  generic  ASCM  accomplishes 
its  search  in  a  matter  of  seconds,  making  search  analogous 
to  taking  a  snapshot  photograph  ( an  aggregated  detection 
process).  Consequently,  target  motion  is  treated  as  negli- 
gible throughout  successive  searches  by  individual  ASCM;  a 
reasonable  assumption  if  the  interarrival  times  between 
missiles  within  the  salvo  are  small.  The  assumption  is  made 
that  each  of  the  reliable  ASCM  in  the  salvo  detect  the  same 
set  of  targets.  This  is  because  of  their  identical  search 
programming,  the  Definite  Range  Law  of  Detection,  and  a 
short  firing  interval. 

Once  a  target  or  number  of  targets  have  been 
detected,  the  missile  enters  an  acquisition  phase  in  which  a 
decision  to  acquire  a  particular  target  is  made.    The  model 
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associates  a  conditional  probability  of  acquisition,  an 
input,  with  each  target  ( HVT  and  one  or  more  LVTs)  detected 
by  the  ASCM.  Implicit  in  the  detection  and  acquisition 
modeling  is  the  assumption  that,  for  purposes  of  acquisi- 
tion, detected  targets  are  processed  independently  and  in 
sequence  from  near  to  far  within  the  missile  seeker  sensor's 
range  gate.  Therefore,  acquisition  of  a  particular  target 
is  influenced  by  two  factors  following  detection:  (a)  its 
relative  position  within  the  range  gate,  and  (b)  its  associ- 
ative conditional  probability  of  acquisition.  An  ASCM 
detection/acquisition  system  with  a  multi-target  discrimina- 
tion ability  has  a  higher  probability  of  acquiring  the  HVT 
than  that  of  acquiring  a  LVT  nearby.  The  model  represents 
this  discrimination  by  incorporating  multiple  Sweep  Widths 
and  conditional  acquisition  probabilities  associated  with 
the  ASCM  system  and  target  set.  In  practice.  Sweep  Widths 
and  conditional  acquisition  probabilities  may  be  difficult 
to  estimate.  Possible  sources  are  ASCM  system  simulations, 
observational  test  firings,  and  expert  opinion. 

One  final  assumption  is  made  concerning  the  nature 
of  the  offensive  firing.  Destruction  of  defensive  firepower 
is  the  initial  result  of  missile  hits  on  the  HVT  and  her 
escorts.  Target  sinkings  occur  long  after  the  ASCM  salvo 
has  run  its  course.  Hence,  although  hits  on  secondary 
targets  degrade  defensive  firepower,  total  target  elimina- 
tion is  not  likely  to  occur  during  the  search  and  detection/ 
acquisition  phase  of  the  ASCM  attack  ( ie. ,  an  HVT  or  LVT 
that  is  present  for  detection  by  the  first  missile  in  the 
salvo  will  be  present  for  the  last  missile  to  detect). 
3.   Modeling  the  Target  Formation 

The  target  formation  is  modeled  as  a  single  HVT 
within  a  bounded,  spatially  homogeneous  Poisson  random  field 
of  escorts.  The  Poisson  field  completely  characterizes  the 
number  of  LVT  in  the  field   ( a  Poisson  random  variable  given 
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by  the  field  density  parameter)  and  the  placement  of  LVT 
within  the  field  (uniformly  randomly  distributed  conditioned 
on  the  number).  It  may  seem  that  the  Poisson  field  is 
better  suited  for  modeling  a  scattering  of  merchant  shipping 
in  a  sea  lane,  rather  than  a  geometrically  rigid  Anti-air 
Warfare  (AAW)  screen.  In  reality,  screens  are  not  rigid: 
escort  ships  are  randomly  patrolling  sectors,  exchanging 
stations,  detaching  to  form  Anti-Submarine  Warfare  (ASW) 
groups,  etc.  ,  making  for  a  dynamic  structure  in  contrast 
with  the  mental  image  of  a  rigidly  symmetric  AAW  screen.  A 
desirable  feature  of  the  Poisson  field  is  the  notion  that 
the  number  of  LVT  in  the  field  is  a  random  variable.  This 
is  likely  to  be  true  if  the  strike  planner's  intelligence 
estimate  of  enemy  force  composition  is  uncertain  and 
targeting  of  secondary  LVT  is  unavailable.  The  average 
number  of  escorts  accompanying  the  HVT  and  the  area  of  the 
Poisson  field  can  be  estimated  with  the  aid  of  intelligence 
reports,  electronic  support  measures,  and  subjective  prob- 
ability encoding,  thereby  furnishing  a  reasonably  good  esti- 
mate of  the  field  density  parameter.  For  example,  a  field 
of  8  escorts  uniformly  distributed  within  a  circle  of  radius 
12  nautical  miles  about  the  HVT  has  a  density  of  0.0177  LVT 
per  square  nautical  mile.  As  will  be  seen,  such  a  figure 
can  be  utilized  to  select  a  nearly  optimum  range  gate 
setting. 

B.   ANALYTICAL  DEVELOPMENT 
1.   Geometry 

A  few  assumptions  concerning  the  orientation  of  the 
Poisson  field  of  escorts  and  the  effective  search  area  of 
the  missiles  are  required  in  order  to  obtain  analytical 
results: 

•  The  HVT  is  centrally  located  within  the  escort  field. 

•  The  escort  field  area  is  large  compared  to  the  ASCM 
search  area  (  ie.  ,  the  search  area  is  wholly  contained  in 
the  field). 
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The  model's  analytics  reference  a  standard 
xy-coordinate  system  oriented  so  that  the  y-axis  lies  on  the 
missile  flight  path.  The  origin  represents  both  the  esti- 
mated mean  of  the  HVT  density  in  residence  above  the 
xy-plane  and  the  center  of  the  Poisson  field  of  escorts. 
The  effective  search  area  of  the  missiles  is  dimensioned  by 
range  gate  (measured  along  the  y-axis)  and  Sweep  Width 
(measured  along  the  x-axis).  It  is  important  to  remember 
that  the  effective  search  area  can  be  a  group  of  rectangles 
with  different  Sweep  Widths  for  different  target  classes 
when  a  distinction  is  made,  in  which  case  each  rectangle  of 
effective  search  has  a  common  y-dimension  (range  gate)  and  a 
unique  x-dimension  (Sweep  Width).  Figure  2.1  is  a  sketch  of 
model  geometry. 

2.   Single  Shot  Penetrator 

For  illustrative  simplicity,  the  HVT  true  location 
density  at  the  time  of  missile  search  is  assumed  to  be 
circular  normal.  Appendix  A  contains  a  time-dependent 
general  bivariate  normal  HVT  location  density  as  an  alterna- 
tive. 

Let  the  random  variable  Y  represent  the  y-coordinate 
of  the  HVT's  true  location,  measured  from  its  estimated 
location  ( ie.  ,  the  origin).  In  the  language  of  Bayesian 
statistics,  the  distribution  of  Y  is  the  posterior  distribu- 
tion of  the  HVT's  true  y-coordinate,  given  a  point  estimate 
of  its  y-coordinate  of  y  =  0.  In  the  following  discussion, 
y  refers  to  a  possible  value  for  the  random  variable  Y. 
(The  same  interpretation  holds  for  the  random  variable  X, 
the  HVT's  true  x-coordinate. ) 

The  single  shot  calculation  requires  the  following 
input  constants  and  probabilities: 

A  =  LVT  escort  field  density. 
W  =  Sweep  Width  governing  LVT. 
W,  =  Sweep  Width  governing  HVT. 
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Figure  2.  1    MODEL  GEOMETRY. 
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d^  =  Near  edge  of  range  gate  (y  coordinate). 

d2  =  Far  edge  of  range  gate  (y  coordinate). 

r  =  Missile  reliability. 

p   =  Conditional   probability  of   acquiring  a 

LVT  given  LVT  is  detected. 

P   =  Conditional  probability  of  acquiring  the 
h 

HVT  given  HVT  is  detected. 

Assume  that  the  HVT  is  contained  in  the  ASCM's 
effective  search  area.  Since  detected  targets  are  processed 
from  near  to  far  within  the  missile  seeker  sensor's  range 
gate,  HVT  acquisition  requires  that  no  LVT  has  been  acquired 
( ie.  ,  selected  and  homed-on)  in  the  region  W  (y-d^ )  .  By 
virtue  of  the  Poisson  field  assumption  this  is  equivalent  to 
experiencing  a  free  path  of  areaVJ  Cy-di)  in  a  in  a  two- 
dimensional  filtered  Poisson  process  of  rate  XP  . 

Define  the  following  events: 

D:   HVT  acquisition. 

E:   Filtered  escort  free  path  to  the  HVT. 

The  probability  of  HVT  acquisition  is  expressed  as 

d2  W  /2 
PCD)  r  /   /^      rP   P(E|Y=y)  f(x,y)  dxdy,       (2.1) 

where 

-(x2+y2)/2a2 
f (x,y)  =  e  1   ,  (2.2) 


2  TT  a 


2 


is  the  probability  density  function  describir 3  the  location 

of  the  HVT,  and 

-XW  P  (y-di) 
P(E|Y=y)  =  e    ^  ^       ,  (2.3) 

is  the  conditional  probability  of  event  E,  given  the  true 
range  coordinate  of  the  HVT.   In  terms  of  these  components, 
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d2  W  /2         _xw  P  (y-di)   -(x2+y2)/2a2 
PCD)  =  /   /^      rP   e    ^  ^        e  dxdy , 

di  -W,  /2     ^  .       2 


XW  P  di 

rP^  e   ^^   {$(W^/2a)-  $(-W^/2a)} 
h  h  h 


•2  /o^2 


d2      -^^oP^y    -y^/2a' 

X  /    e         e         dy 


/2  7Ta 


The  integral  is  evaluated  by  completing  the  square  in  the 
exponent  and  adjusting  the  limits  of  integration. 

Let  K  =  XW  P„,  z  =  y±ii^     ^    ^2  =  ^; 


Kd^  +'2(  Ko) 


2 


PCD)  =  rP,  e           A$  {^Cdo+^a^)  -  $  Cd  i +Ka^  )  }, 
h  X    — ^ — ■'■ 


a  a 


1^  r  ^  ^  ^  2 


Kd^  +i$C  KO) 

=  rP,  e  A$A$',  C2.U) 

h  X   y 


where 


b 
A$  =  /   ct)Cu)  du 
a 


is  the  integral  of  the  standard  normal  probability  density 
function  evaluated  at  the  upper  and  lower  limits  shown. 
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3.   Multiple  Penetrators 

The  probability  that  at  least  one  of  m  penetrating 
missiles  acquires  the  HVT  is  computed  by  conditioning  on  the 
number  of  LVT  in  W.(y-di)  of  the  effective  search  area.  The 
solution  is  reached  by  integrating  the  HVT  density  over  the 
effective  search  area,  after  first  removing  the  condition  on 
the  number  of  LVT  "in  the  way". 

Let  N,  a  random  variable,  be  the  number  of  missiles 
that  acquire  the  HVT.  The  conditional  probability  that  at 
least  one  of  m  missiles  acquires  the  HVT  is: 


1  ^m 


P(N>0|L=i,  Y=y)  =  1  -  (1  -  rP(l-P  )  ) 


(2.5) 


where  L  is  the  number  of  LVT  in  W  (y-di). 


The  condition  on  LVT  is  removed  using  the  knowledge  that  the 
assumed  distribution  of  LVT  escorts  is  Poisson  with  spatial 
density  parameter  X  .  Expansion  of  the  binomial  term  and 
rearrangement  shows  that 


P(N>0 lY^y)  =  1  - 


;  [ 

i  =  0  L 


(1  -  rP^d-P)"-)"" 
h     £ 


X  e 


-^W^^y-^l^  (XWCy-di))^ 


«>    m 

1  -   Z    E 
i=0  k=0 


\l) 


(-l)^(rP,  (1-P  )^)^ 
n     £ 


X  e 


-^^/y-^l^  (AW^(y-di))^ 
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m 
1  -   Z 
k=0 


\k/  (-rP,  )   e 

h 


i=0        ^  i 


T  I 


m 
1  -   E 
k=0 


\k/  (-rPi_)   e 


(1-P^)^ 


(2.6) 


Finally,  to  remove  the  condition  on  the  location  of  the  HVT, 
equation  (2.6)  is  multiplied  by  the  HVT  location  density  and 
integrated  using  the  same  technique  employed  for  the  single 
shot  calculation. 


Let  5(k)  =  XW^d  -  (1-P^)^); 


d2  W  /2  m 

P(N>0)  =  1  -  /  /  I 

di  -W^/2  k=0  I 

^         n 


(k) 


(-rP,  )^  e 

h 


-5(k)(y-di) 


-(x2+y2)/2a2 
X  e  dxdy 

2TTa2j 


=  1  -  A' 


X 


m 

E 

k=0 


il)  c- 


6(k)d 


rP^)   e 

h 


d2    -6(k)y   -y2/2a2 
X  /    e        e         dy 

^1  /27a 
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m 


6(.k)(li+h(    6(k)a    ) 


1    -    A$         S  Ik/     (-rP,  )       e  A*,.'  5 

^    k=0  ^  ^ 


m 
k  =  l 


il) 


2 


6(k)di+i^(    6(k)a    ) 

A$         E  ik^     (-r^P,.)       e  A<I>    '  , 

X   ,     ^  h  y 


(2.7) 


where 

A$    '     =    $(do+6(k)a2)    -    $ ( dj+6(k)o2 ) . 
V  — 7^ 

Expressions  (2.4)  and  (2,7)  allow  the  analyst  to  determine 
an  optimal  near  edge  of  range  gate  setting,  di ,  ie.  one  that 
maximizes  HVT  acquisition  probability,  given  the  various 
operational  parameters  of  the  ASCM. 

4.   The  Fall  of  Multiple  Penetrators 

More  information  can  be  obtained  from  the  condi- 
tional distribution  of  the  fall  (distribution)  of  multiple 
penetrators,  given  that  the  HVT  is  contained  in  the 
missile's  effective  search  area.  Fall  of  shot  terminology 
(e.g.,  dispersion)  is  normally  associated  with  the  firing 
assessment  of  ballistic  projectiles,  but  can  be  adopted  for 
ASCM  if  one  can  estimate  the  probability  of  "absorption"  of 
an  ASCM  into  each  target  class.  Simply  stated,  a  missile 
which  penetrates  the  HVT  area-defense  will  either  have  the 
HVT  in  its  projected  search  area  or  not.  If  not,  it  is  no 
longer  a  threat  to  the  HVT,  though  there  is  a  probability, 
usually  small,  that  it  will  home  on  an  escort.  For  the 
missiles  that  are  a  threat  to  the  HVT,  it  is  possible  to 
estimate  how  they  are  distributed  between  HVT,  escorts,  and 
misses. 
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The  acquisition  decision  of  a  single,  conditionally 
independent  penetrator  is  likened  to  a  single  trial  of  a 
multinomial  experiment  which,  given  the  HVT  is  in  the  effec- 
tive search  area,  has  as  its  possible  outcomes: 


ASCM  acquires  a  LVT  on  the  near-side  of  the  HVT. 
ASCM  acquires  the  HVT. 

ASCM  acquires  a  LVT  on  the  far-side  of  the  HVT. 
ASCM  fails  to  acquire  a  target. 


The  probability  formulas  of  the  above  events  are 
initially  stated  with  the  condition  of  the  HVT's 
y-coordinate  in  the  effective  search  area. 

Let  Ad  =  d2  -  di,  and  condition  on  Y=y;  then 

-K (y-d  Y  ) 
P(A|Y=y)  =    r(l  -  e         ) ,  (2.8) 

-K (y-d J  ) 
P(B|Y=y)  =  rP^  e         ,  (2.9) 

-K (y-dj  )        -K(d2-y ) 
P(ClY  =  y)  =  r(l-P^)(e         )(1  -  e        .)  ,    (2.10) 

-KAd 
P(D|Y  =  y)  =  (1-r)  +  r(l-P,  )  e     .  (2.11) 

Derivations  of  the  mean,  variance,  and  covariance  of  the 
multinomial  distribution's  binomially  distributed  marginal 
distributions  are  given  in  many  basic  probability  texts. 
Removing  the  condition  of  the  y-coordinate  of  the  HVT  within 
the  effective  search  area  and  normalizing  with  respect  to  y 
produces  the  following  marginal  expectations  of  the  condi- 
tional multinomial  events  ( ie. ,  the  expected  number  of 
missiles  absorbed  by  each  target  class,  given  the  HVT  is  in 
the  effective  search  area). 
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Let  M.  be  the  random  number  of  penetrators  that  are 
"absorbed"  into  cell  i,  where  M .  +  M_  +  M„  +  M^  =  m .  The  ex- 
pected  value  of  M.  is  mp . ,  where 


p  =    P(A)  =  r  -  (1/A<|)  )  /    re 
^  y   di 


-K(y-di)   -y2/2a' 


dy 


2-n  o 


Kd I +^( Ka ) ^ 

r(l  -  e  )A$  '  /A$ 

y   y 


(2.12) 


p^  =  P(B)  =    (1/A$  )  /    rP,  e 
^  y   di     ^ 


-K(y-di)   -y2/2a2 

e         dy 


27Ta 


=  rP,  e  A$  '/A$  ; 

h  y    y 


(2.13) 


p   z  p(C)  =  (1/A<D  )  / 

^  y  di 


r(l-P^) (e 

h 


-K (y-d^ )     -KAd 
-  e     ) 


-y2/2a- 


X  e 


dy 


2-n  o 


K:di+ig(Ka)2 
r(l-P,)(e  A$  '/A$ 

h  y    y 


-KAd 
e    ); 


(2.14) 


d  9  r 
p^  =  P(D)  =  (1-r)  +  (1/A$  )  /' 

D  y  ^ 

^   di  L 


-KAd 


r(l-P^)  e 
h 


-y2/2a2 
X  e         dy 


/2  7Ta 
-KAd 


(1-r)  +  r(l-P,  )  e 
h 


(2.15) 
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The  general  procedure  for  computing  a  marginal  variance  is 
shown  below,  followed  by  the  equation  of  the  variance  of 
event  B;  the  number  of  HVT  acquisitions,  given  the  HVT  lies 
within  the  effective  search  area.   In  general, 

Var(Y)  =  E{  Var(Y|X)  }  +  Var{  E(Y1x)  }. 
Therefore , 

Var(M.)  =    E{mp.(l  -  p.(y))}  +  Var{mp.(y)}, 


=  E{mp. (1  -  p.  (y)) }  +  m2 


fEfp? 


(y)}  -  (E{mp.(y)}) 


=  m(m  -  l)E{p?(y)}  +  mE{p.(y)}  -  m2(E{p.(y)})  , 


m(m  -  l)E{p?(y)}  +  mp.  -  (mp.)  . 


(2.16) 


Equations  (2.8)  through  (2.9)  can  be  substituted  into  equation 
(2.16)  in  place  of  p.(y)  and  the  corresponding  expressions 
(2.12)  through  (2.15)  substituted  in  place  of  p . .   The  expect- 
ation of  p?(y)  remains  to  be  calculated  and  can  be  accomplished 
using 

d2  /  d2 

E{g(p.(y))}  =  /    g(p.(y))  (})(y)  dy/  /    ^(y)  dy 
di      ^  /   di 


the  expectation  of  a  function  of  p.(y),  given  the  HVT  is  in 
the  missile's  effective  search  area  (hence  the  normalization 
quotient)  . 

Calculation  of  the  variance  of  M„ ,  the  number  of 

D 

ASCM  that  acquire  the  HVT,  given  the  HVT  is  contained  in  the 
missile's  effective  search  area  is  as  follows; 
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Var(M^)    =   m(m   - 


d2 

1)(1/A$    )    / 


h 


-2K(y-di) 


-y2/2a^ 
X    e  dy 


/2 


Tia 


mp. 


m(m   -    Dr^p    2    e 
n 


2(Kdi+(Ka)2) 


(mp^)     , 


y  y 


mp. 


(mp^)     , 


(2.17) 


where 


A$ 


y 


(d:>  +  2Ka^) 


(di +2Ka^) 


It  is  important  that  the  effect  of  the  range  gate 
decision  variable  on  HVT  acquisition  probability  can  be 
explored  using  the  model's  analytical  equations.  An  illus- 
trative numerical  example  is  presented  in  the  next  chapter. 
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III.  OPTIMAL  RANGE  GATING 

The  near  edge  of  the  missile  seeker  sensor's  range  gate, 
d^ ,  is  an  influential  decision  variable.  As  the  range  gate 
increases  in  length,  more  targets  are  likely  to  be  detected 
and  processed  for  acquisition,  making  the  ASCM's  task  of 
selecting  the  HVT  more  difficult.  As  the  range  gate 
narrows,  the  probability  of  capturing  the  HVT  in  the  effec- 
tive search  area  (  ie.  ,  detecting  the  HVT)  decreases.  If  the 
far  edge  of  the  range  gate,  d2  ,  has  been  chosen  such  that 
the  subsequent  choice  of  d^  brackets  the  mean  of  the  HVT 
targeting  density  along  the  missile  flight  path,  equation 
(2.7)  can  be  used  to  solve  for  a  value  of  d^  that  maximizes 
the  probability  that  at  least  one  penetrating  missile 
acquires  the  HVT,  given  m  penetrators  are  realized.  If  an 
estimated  distribution  of  the  number  of  penetrators  were 
available,  the  law  of  total  probability  could  be  used  to 
obtain  an  unconditionally  optimal  or  "weighted"  value  of  d^ 
based  on  the  conditionally  optimal  values  of  d^  obtained 
from  graphing  equation  (2.7)  for  m  penetrators,  m  = 
l,2,...,s,  where  s  is  the  total  salvo  size.  Or,  in  the 
situation  where  only  an  expected  number  of  penetrators  is 
available,  equation  (2.7)  can  be  graphed  once  for  the 
optimal  value  of  d,  based  on  the  expected  value  of  m.  The 
numerical  example  that  follows  illustrates  how  the  optimal 
choice  of  d  is  affected  by  ASCM  multi-target  discrimination 
ability,  the  number  of  penetrators,  and  HVT  targeting  accu- 
racy. 

The  example  uses  the  following  hypothetical  input  param- 
eter values: 

^  =  3.0  Nm 

X    =   0.  0177  LVT/Nm2 


W,  =  6  Nm 
n 
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Wj^  =  6  Nm 

Py^  =  0.97 

P,  =0.4 

^  =  0.  99 

d2  =  +5. 0  Nm 


Figure  3. 1  is  a  plot  of  a  single  shot  penetrator's  HVT 
acquisition  probability  versus  di ,  where  zero  on  the  hori- 
zontal axis  is  the  expected  location  of  the  HVT  along  the 
flight  path.  The  range  gate  [  d^  , d2 ]  =  [-5.1, +5.0]  maximizes 
the  single  shot  probability  of  HVT  acquisition.  Figure  3.2 
contrasts  the  curve  of  Figure  3.1  (curve  A)  with  a  curve 
generated  from  identical  input  parameter  values  excepting 
conditional  acquisition  probabilities,  which  have  been  set 
equal  to  0.97  (curve  B).  The  relatively  gradual  ramp  of 
curve  A  is  representative  of  a  "more  forgiving"  missile 
( ie.  ,  one  that  has  a  better  multi-target  discrimination 
ability) . 

Figure  3. 3  demonstrates  the  increase  in  HVT  acquisition 
probability  from  firing  a  salvo  and  realizing  four  penetra- 
tors  versus  firing  a  single  shot  with  the  hope  of  realizing 
one.  While  it  is  difficult  to  spot  a  noticeable  difference 
in  d  associated  with  the  maxima  of  each  curve,  the  actual 
difference  is  about  1. 7  nautical  mile,  with  the  four  pene- 
trators  having  the  "longer"  optimal  range  gate.  The  result 
is  a  consequence  of  maximizing  the  probability  that  at  least 
one  missile  acquires  the  HVT.  The  marginal  "single  shot" 
probability  of  HVT  acquisition  of  each  missile  is  effec- 
tively reduced  while  the  chance  at  least  one  acquires  the 
HVT  is  increased.  This  phenomenon  becomes  more  apparent  in 
Figure  3. 4,  ten  penetrators  versus  a  single  shot  penetrator. 

Figure  3.5  is  a  plot  of  number  of  penetrators  versus  HVT 
acquisition  probability  where  the  conditionally  optimal 
range  gate  has   been  used  at  each  increment   of  the  indepen- 
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dent  variable.  The  marginal  return  of  the  first  few  pene- 
trators  is  striking.  Figures  3.7  and  3.8  demonstrate  the 
marked  effect  that  an  increasing  HVT  location  error  sigma 
has  on  HVT  acquisition  probability  in  both  the  single  shot 
penetrator  and  salvo/multiple  penetrator  cases.  As  sigma 
increases,  it  becomes  apparent  that  the  simple  firing 
strategy  may  very  well  be  out-performed  by  a  strategy  that 
involves  firing  a  spread  of  missile  salvos  (assets  permit- 
ting), an  interesting  hypothesis  for  future  research. 

Table  I  summarizes  conditional  fall  of  shot  expectation 
and  standard  deviation  for  a  single  shot  penetrator  and 
various  numbers  of  multiple  penetrators  over  selected  values 
of  sigma.  Again,  the  fall  of  shot  for  events  A,  B,  C,  and  D 
is  conditioned  on  the  HVT's  containment  and  location  within 
the  missile's  effective  search  area.  As  sigma  increases, 
the  optimal  range  gate  opens,  the  probability  of  containing 
the  HVT  in  the  effective  search  area  decreases,  and,  in  the 
case  of  salvo/multiple  penetrators,  the  conditional  fall  of 
shot  expectation  becomes  more  evenly  distributed  between  the 
HVT  and  the  LVT  on  the  near  -  side.  It  should  be  noted  that 
the  total  expected  threat  to  the  escorts  cannot  be  deter- 
mined by  simply  adding  the  expectation  of  events  A  and  C, 
because  there  is  a  chance  that  some  escorts  may  lie  in  the 
missile's  effective  search  area  though  the  HVT  does  not. 
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Figure   3.1         SINGLE    SHOT   P(HVT   ACQ.  )    VS.     RANGE   GATE. 
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Figure  3.2    THE  EFFECT  OF  MULTI-TARGET  DISCRIMINATION. 
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Figure   3.3        ONE   VS.    MANY. 
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Figure    3.4        RELAXING  THE   RANGE   GATE   DECISION  VARIABLE. 
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Figure  3.5    INCREASING  PENETRATORS, 
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Figure    3.6         TARGETING  ACCURACY  AND   THE    SINGLE    SHOT. 
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Figure  3.7    TARGETING  ACCURACY  AND  MULTIPLE  PENETRATORS. 
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TABLE  I 
CONDITIONAL  FALL  OF  SHOT 

SINGLE  SHOT  PENETRATOR 

Expectation  ( Std,  dev.  )  of  events; 
a  "  (dj  ,d2  )    P(HVT)*      A        B        C        D 

1.5     -3.0, +5.0    0.932     0.119    0.844    0.004    0.031 

(0.  362) 

3.0  -5.1, +5.0         0.619  0.188         0.777         0.004        0.029 

(0.  415) 

9.0         -12.1, +5.0         0.162  0.315         0.654        0.005         0.024 

(0.  475) 


FOUR  PENETRATORS 

Expectation  ( Std.  dev. )  of  events: 
a  '  (d^  ,d2  )    P(HVT)*      A       B        C       D 

1.5     -4.0, +5.0    0.950     0.614    3.24     0.019    0.121 

(0.  801) 

3.0  -6.  8, +5.0         0.642  0.947         2.92  0.018         0.111 

(  1.  000) 

9.0         -16.0, +5.0         0.175  1.57  2.31  0.022         0.088 

( 1. 310) 


TEN  PENETRATORS 

Expectation  ( Std.  dev. )  of  events: 
a  {d^,d^)         P(HVT)*      A       B       C       D 

1.5     -4.  8, +5.0    0.953     1.81     7.84     0.046    0.295 

(1.38) 

3.0  -8.  3, +5.0         0.648  2.81  6.87  0.043         0.268 

(1.89) 

9.0         -19.1, +5.0         0.181  4.54  5.19  0.053         0.206 

(2. 66) 

*  P(HVT)  =  "The  prob.  of  HVT  containment  in  the  search  area". 
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IV.  SUMMARY 

The  number  of  ASCM  area  defense  penetrators  necessary  to 
achieve  a  required  minimum  number  of  HVT  acquisitions  is  a 
random  variable.  Tables  or  nomographs  containing  number  of 
penetrators,  optimal  range  gate,  expected  number  of  HVT 
acquisitions,  and  standard  deviation  of  HVT  acquisitions 
versus  escort  field  density  could  be  generated  for  a  set  of 
different  posterior  distributions  describing  HVT  location 
over  a  range  of  Poisson  field  densities.  It  would  be  neces- 
sary to  accompany  these  "conditional"  estimates  with  a  prob- 
ability statement  of  HVT  containment  in  the  effective  search 
area.  The  model  in  its  present  form  serves  best  as  an 
example  of  how  a  probability  model  incorporating  a  rela- 
tively simple  analytic  form  can  be  used  for  OTH  strike  plan- 
ning. Complex  firing  situations  and  target  coverage 
problems  involving  multiple  launches  from  a  consortium  of 
strike  groups  could  be  addressed  by  expanding  the  simple 
analytical  framework  of  the  model.  The  equations  of  a  well- 
tailored  ( and  well-validated)  follow-on  model  could  be 
contained  in  an  interactive  computer  program  which  accepts 
periodic  targeting  updates  and  in  turn  furnishes  the  strike 
planner  with  new  estimates  of  salvo  size  and  optimal  range 
gating  to  achieve  a  desired  probability  of  HVT  acquisition 
based  on  the  number  of  area  defense  penetrators. 
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APPE^DIX  A 
UPDATING  THE  HVT  LOCATION  DENSITY 

The  time  delay  measured  from  the  generation  of  the  HVT 
location  density  to  the  time  the  missile  commences  its 
search  will  often  allow  considerable  movement  of  the  enemy 
formation.  A  postulated  bivariate  normal  target  velocity- 
time  vector  can  be  used  to  revise  the  "old"  HVT  density. 
The  approach  does  not  win  "something  for  nothing"  since  an 
assumption  concerning  enemy  intentions  is  required. 
Constant  target  motion  is  a  reasonable  assumption  if  (a)  the 
enemy  has  not  been  alerted  by  the  opposition' s  tracking/ 
targeting  methods  and  (b)  timely  early-warning  of  the 
impending  strike  is  unlikely.  A  Brownian  motion  or  fleeing 
datum  assumption  would  be  more  appropriate  in  the  alerted 
enemy  situation.  Dispersion  of  the  enemy  formation  is  also 
a  possibility  in  the  event  of  a  probable  nuclear  exchange. 
The  constant  motion  model  developed  below  may  prove  valid 
when  offensive  targeting  is  reasonably  passive  and  covert. 

Let  OX,  the  HVT's  vector  of  movement  to  its  future  posi- 
tion measured  from  the  mean  of  the  "old"  HVT  location 
density,  be  the  resultant  of  the  sum  of  two  independent 
bivariate  normal  vectors:  a  position  error,  OX' ,  and  a 
velocity-time  vector  X'X,  (see  Figure  A.  1 ) . 
The  distribution  of  OX'  is  given  by  the  HVT  density, 

OX'^^BYNCy^,  Py'  ^x  '  ^Y '  ^XY"^ 
while  X'X  is  distributed  as  follows, 

X'X '^-BVN(pyt,  y^t,  a^t^,  a2t2,  p^^) 
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ox  =    OX'  +  X'X 


Figure  A.  1    THE  SUM  OF  TWO  INDEPENDENT  RANDOM  VECTORS. 
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The  distribution  of  OX  is  gleaned  from  its  moment  gener- 
ating function  formed  by  the  product  of  the  individual 
moment  generating  functions  of  OX'  and  X'X.  The  moment 
generating  function  of  the  bivariate  normal  distribution  is: 


X,Y 


m,,  ,,(si,S2)  =  exp{siy^  +  S2yY  ^    hislc^   +  2psiS2a^a 

+  s|a2)}. 

Therefore, 

1.  r  c-  2  ^  2 


Y 


m 


+  s^a^) } , 
m,,,y(si,S2)  =  expCs^y  t  +  S2y„t  +  h{s\o^t^ 


+  2p^^SiS2ayta^t  +  s^a2t2)}, 


m^,,(si,S2)  =  m^„,(si,S2)  m  ,„(si,S2) 


and, 

1q^v=,1  ,^2^  -  -OX'  '"^  '"^'  '"X'X 

=  exp{si(y^  +  y^t)  +  S2(yY  +  y^t^ 

•   -      +   m    sf(a2  +  a^t^)  +  2SlS2(p^^Y^X^Y 

By  inspection, 

OX'^BVNCy^  +  y^t,  y^  +  ^v^'  °X  "^  "^U^^'  '^Y  "^  '^V^^' 

PXY^X^Y  ""  PUVV^V^^' 
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The  distribution  of  OX  can  be  used  to  adjust  the  launch 
bearing  of  the  ASCM  salvo  and  update  HVT  postional  uncer- 
tainty to  the  estimated  time  of  missile  search.  Model  geom- 
etry becomes  more  complicated  when  a  general  bivariate 
normal  HVT  location  density  is  incorporated  instead  of  a 
circular  normal  one  since  a  unique  orientation  of  the  effec- 
tive missile  search  area  results. 
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